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Abstract 

Rescaled Mellin-type transforms of the exponential functional of the Bourgade-Kuan-Rodgers 
statistic of Riemann zeroes are conjecturally related to the distribution of the total mass of the 
limit lognormal stochastic measure of Mandelbrot-Bacry-Muzy. The conjecture implies that a 
non-trivial, log-infinitely divisible probability distribution is associated with Riemann zeroes. 

For application, integral moments, covariance structure, multiscaling spectrum, and asymptotics 
associated with the exponential functional are computed in closed form using the known mero- 
morphic extension of the Selberg integral. 
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1 Introduction 

In this paper we contribute to the literature on the statistical distribution of Riemann zeroes in the 
mesoscopic regime. The study of the values of the Riemann zeta function in the mesoscopic regime 
was pioneered by Selberg ifSTI . |[5^ and then extended to the distribution of zeroes by Fujii 1231, 
Hughes and Rudnick l32l . Bourgade fT4ll . and most recently by Bourgade and Kuan ITSl . Rodgers 
l53l . and Kargin l35l . Assuming the Riemann hypothesis (except for Selberg), these authors rigor¬ 
ously established various central limit theorems for the distribution of Riemann zeroes. The principal 
technical tools that were used to obtain these theorems were Selberg’s formula for explicit 
formulas of Guinand and Weil, and certain moment calculations. Alternatively, beginning with the 
seminal work of Montgomery |@0l, a great deal of progress has been made in formulating precise 
conjectures about the statistical distribution of the zeroes. These conjectures are all motivated by the 
empirical fact that the statistical properties of the zeroes are very close to those of eigenvalues of large 
Hermitian matrices with independent entries, i.e. the so-called GUE matrices, up to small arithmetic 
corrections, and calculations are typically justified by means of semi-classical methods for quantum 
chaotic systems, Keating-Snaith philosophy of modeling the value distribution of the Riemann zeta 
function on the critical line by the characteristic polynomials of certain large random matrices, and 
conjectural forms of the approximate functional equation for the zeta function. Eor example. Berry 
uni calculated the GUE term and arithmetic corrections for the number variance, Bogomolny and 
Keating ifT^ did the same for the pair correlation function, which was later extended to multiple- 
point correlations by Conrey and Snaith ifTTl and Bogomolny and Keating |[T3l . Keating and Snaith 
Il3^ calculated the moments of the zeta function on the critical line, Conrey et. al. ifT^ formulated the 
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ratios conjecture for its average values, Farmer et. al. iBOll and Fyodorov and Keating estimated 
the magnitude of its extreme values on the critical line, to name a few. 

In this paper we conjecture a mod-gaussian limit theorem associated with the distribution of Rie- 
mann zeroes in the mesoscopic regime by combining the approach of Bourgade and Kuan with our 
previous work on the limit lognormal stochastic measure (also known as lognormal multiplicative 
chaos) and Selberg integral. Bourgade and Kuan and Rodgers independently proved that a class of 
linear statistics of Riemann zeroes converge to gaussian vectors and, most importantly, computed 
the covariance of the limiting vector explicitly. The starting point of our approach is that this lim¬ 
iting gaussian vector approximates the centered gaussian free field when the statistic is based on a 
smoothed indicator function of subintervals of the unit interval, and it approximates the gaussian free 
field plus an independent gaussian random variable when the statistic is based on a smoothed indi¬ 
cator function of subintervals of a particular unbounded interval. The limit lognormal measure is 
defined as a limit of the exponential functional of the gaussian free field, hence by taking the Mellin 
transforrr0 of the exponential functional of the Bourgade-Kuan-Rodgers statistic in an appropriately 
rescaled limit, we conjecturally obtain the Mellin transform of the total mass of the measure. In a se¬ 
ries of papers, cf. 1431, 1441 . ll45]|. 1481 . we investigated the total mass and made a precise conjecture 
about its probability distributiono The positive integral moments of the total mass are known to be 
given by the classical Selberg integral. In 1451 we rigorously constructed a probability distribution 
(called the Selberg integral distribution) whose nth moment coincides with the Selberg integral of 
dimension n and conjectured that this distribution is the same as the distribution of the total mass. 
Thus, the main result of this paper is a conjecture that particular rescaled limits of two Mellin-type 
transforms of the exponential functional of the Bourgade-Kuan-Rodgers statistic corresponding to 
a smoothed indicator function of certain bounded or unbounded intervals coincide with the Mellin 
transform of the Selberg integral distribution. In BSl . Il49]l . |[50l we established many properties of 
the Mellin transform so that we can make a number of precise statements about our rescaled limit as 
corollaries of the main conjecture. The type of rescaling and convergence that we use in this paper 
is closely related to the rescaling that was used by Keating and Snaith ll^ and Nikeghbali and Yor 
B^ . formalized by Jacod et. al. If^ in their theory of mod-gaussian convergence, and significantly 
extended in recent publications of Feray et. al. |[^ and Meliot and Nikeghbali 

The main technical innovation of our work is the explicit use of the gaussian free field, limif 
lognormal measure, Selberg integral, and Selberg infegral disfribufion in fhe confexf of fhe sfafisfical 
disfribufion of Riemann zeroes. We note fhaf fhe Selberg infegral and Selberg infegral disfribufion 
previously appeared, respecfively, in conjecfures of Keafing and Snaifh abouf fhe momenfs and 
of Fyodorov and Keafing ll^ abouf exfreme values of fhe Riemann zefa funcfion on fhe crifical line. 
These conjecfures are based on fhe analogy befween fhe value disfribufion of l^{\/2 + it) and fhaf of 
fhe characferisfic polynomials of certain large random mafrices. Our conjecfure deals insfead wifh fhe 
zeroes of fhe zefa funcfion on fhe crifical line and is based on fhe convergence of particular sfafisfics 
of fhe zeroes fo fhe gaussian free field or ifs centered version. In parficular, we prove fhaf our sfafisfics 
or, equivalenfly, particular integrals of Imlog(^(l/2-|-ti) along fhe critical line, exhibif logarifhmic 
correlafions and calculate fhe corresponding covariances explicifly. We believe fhaf fhese calculafions 
are newH 

*It is more natural to define the Mellin transform as JJ” f{x)dx as opposed to the usual JJ” x^^^ fix) dx for our 
purposes. 

^The terms “probahility distribution”, “law”, and “random variable” are used interchangeably in this paper. 

^The idea that Imlog 1/2 -f if) is logarithmically correlated is not new. Keating and Snaith 1361 and later Farmer et. al. 
1^ argued that Imlog (1/2 -f it) can be modeled by the imaginary part of the logarithm of the characteristic polynomial 
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The limit lognormal measure was introduced and reviewed by Mandelbrot |[^ in the context 
of mathematical modeling of intermittent turbulence and constructed explicitly by Bacry and Muzy 
iSl. |4], El, iHTI . Its existence and basic properties follow from the general theory of multiplicative 
chaos of Kahane O. This measure is of significant interest in mathematical physics as it naturally 
appears in a wide spectrum of problems ranging from conformal field theory ifTOl . ll52l and two- 
dimensional quantum gravity IT^ . to statistical mechanics of disordered energy landscapes ll2^ . |l28l, 
im, OOll . to name a few. A periodized version of the limit lognormal measure appears in a random 
energy model 1251 and in the theory of conformal weldings |5|. We also mention a multidimensional 
extension of the measure [H and a recent construction of the critical lognormal multiplicative chaos 
in ifTSll and fS]. One of the most remarkable properties of this measure is that it is stochastically 
self-similar with lognormal multipliers (hence its name), so that the moments of its total mass do not 
scale linearly but rather quadratically, i.e. the measure is multifractal. The aforementioned problems 
in mathematical physics all exhibit multifractal behavior so that the significance of our conjecture 
extends beyond the distribution of Riemann zeroes per se for it suggests that the phenomenon of 
multifractality might have a number theoretic origin in the sense that the distribution of Riemann 
zeroes (conjecturally) provides a natural model for such phenomena. 

We do not have a mathematically rigorous proof of our conjecture and provide instead some exact 
calculations (a ’’physicist’s proof’) that explain how we arrived at it. If one assumes the conjecture 
to be true, the resulting corollaries are mathematically rigorous and their proofs can be found in ll48l 
and Il50l . 

The plan of this paper is as follows. In Section 2 we give a brief review of the key results of 
Bourgade and Kuan and Rodgers and then state our conjecture and its implications. This section does 
not require any knowledge of the limit lognormal measure. In Section 3 we review the limit lognormal 
measure and the Selberg integral distribution. In Section 4 we present a heuristic derivation of our 
conjecture. Conclusions are given in Section 5. 

2 Results 

We begin this section with a brief description of the statistic of Riemann zeroes that was introduced 
by Bourgade and Kuan ifTSl and Rodgers ll5^ (henceforth referred to as the BKR statistic^ following 
the approach and notations of Bourgade and Kuan. The Riemann zeta function is defineqjby 

oo 

C(5) = ^(m-hl)^",Re(5) > 1, (1) 

m=0 

and is continued analytically to the complex plane having a simple pole at 5 = 1. Its non-trivial zeroes 
are known to be located in the critical strip 0 < Re(5) < 1 and, according to the Riemann hypothesis, 
are thought to lie on the critical line Re ( 5 ) = 1/2, cf. Il59l for details. Assuming the Riemann 
hypothesis, we write non-trivial zeroes of the Riemann zeta function in the form {1/2 -|- /y}, 7 G M. 
Let A(t) be a function of t > 0 that satisfies the asymptotic condition 

l<CA(f)<Clogt (2) 

of CUE matrices, and Hughes et. al. ED proved that the latter is logarithmically correlated, thereby conjecturing the 
same about ImlogiJ(l/2 + if). The novelty of our work is the computation of the logarithmic covariance structure of 
Imlog 1/2 -f if) from first principles. 

^We will use the symbol = to mean that the left-hand side is defined by the right-hand side. 
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in the limit t —)• oo^ where the number theoretic notation a{t) <C b{t) means a{t) = o(b{t)Y Let co 
denote a uniform random variable over (1,2), y{t) = X{t){y — (Ot), and define the statistic 

St{f) = £/(7(0) - 2nx\t) I 

given a test function f{x). The class of test functions that was considered in ITSll is primarily 

defined by the condition (/, /) < oo^ where 

(/,§) -Re I\w\f{w)g{w)dw, (4) 

= fi^)s'{y)^og\x-y\dxdy (5) 

plus some mild conditions on the growth of f{x) and its Fourier transform /(w) = I/In J f{x)e~""^ dx 
at infinity and a bounded variation condition (that are satisfied by compactly supported functions, 
by example). We note that St{f) is centereq£l in the limit f —)■ oo as it is well known that the number 
of Riemann zeroes in the interval [f, 2t] is asymptotic to t logt/27r in this limit. The principal results 
of |1T5]| and llssfl and the starting point of our construction are the following theorems. 

Theorem 2.1 (Convergence to a gaussian vector) Given test function fi, ■ ■■ fk in the ran¬ 
dom vector ,St{fk)) converges in law in the limit t ^ to a centered gaussian vector 

(^(/i), • • ■S{fk)) having the covariance 

Coy{S{fi),S{fj)) = {f,fj). (6) 

The second theorem deals with the case of diverging limiting variance. Define 

logf/A(f) 

a,(/)2= j \w\\f{w)\^dw, (7) 

-logr/A(() 

fhen, under fhe assumption that Gt{f) —as t —oo^ 

Theorem 2.2 (Convergence in the case of diverging variance) 

( 8 ) 

where ^(0,1) denotes the standard gaussian random variable with the zero mean and unit variance. 

The significance of fhe condition A(t) <C logt is that the number of zeroes that are visited by / as 
f oo goes to infinity, i.e. Theorems 12. II and 12.21 are mesoscopic central limit theorems. 

The intuitive meaning of the BKR theorems and the statistic St (/) can be established from the 
connection of St{f) with the error term S{t) in the zero counting function N{t). Let N{t) denote the 
number of Riemann zeroes having their imaginary part (“height”) between zero and t. Let the function 
^(f) be defined by 

S{t) = — arg^(l/2 + /t) = — Imlog(^(l/2 + /f). (9) 

71 71 

^Centered means that its expectation is zero. All expectations, covariances, etc in this section are with respect to the 
distribution of £0. 

^Rodgers considered a more restrictive class of test functions and stated the formula for the variance only. 
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Then, the Riemann-von Mangoldt formula states 


N{t) = -(^Imlogr(- + ^) 


^logTT^ +1 + 5(0, 




( 10 ) 

( 11 ) 


in the limit f —> oo. Let u be fixed and X{t) be as in (ID, then the asymptotic in equation (fTTI) implies 
that the number of zeroes in the random interval [(Ot, (Ot + u/X{t)\ satisfies in fhe same limif 


“ “ 2^++ X ^)"+ o[VKt)) • 


X{t), 


2kX (f) 


X{t). 


( 12 ) 


If follows fhaf fhe expected number of zeroes in [cot, (Ot + u/X{t)] is given by fhe leading asympfofic 
Mlogt/27rA(0, whereas S{^(Ot + u/X{t)^ —S{(Ot) gives “fhe error”, i.e. fhe fluclualion of fhe number 
of zeroes in fhis inferval from ifs leading asympfofic. Note fhaf fhe condition in ([D means fhaf fhe 
lengfh of fhe inferval goes fo zero, whereas fhe expected number of zeroes goes fo infinify in fhe limif 
t —oo, i.e. fhe inferval is mesoscopic. Lef Xu{x) denofe fhe indicafor function of fhe interval [0, u], 
fhen fhe corresponding BKR sfafisfic salisfies by ([D and (fT2l) 

5,te,)=A'(«,< + ^)-iV(a»)-«^. (13) 

= 5(mt + ^)-S(m0 + O(l/A(0). (14) 


Hence fhe sfafisfics measures fhe fluclualion of fhe error term over fhe random inferval [(Ot,(Ot + 
u/X{t)]. If is easy fo see from ([71) fhaf fhe corresponding variance has fhe leading asymptotic o^{Xu) ~ 
log(logt/A(f)) /tt^ so fhaf Theorem 12.21 gives us 

Sjc ot + u/X{t)) - 5 ^ ^^ 

^log{logt/X{t))/K 

This special case of Theorem l2.2l is known as Fujii’s cenlral limif Iheorem, cf. |[23l . f2M . If lurns oul 
fhaf fhe inferprefafion of fhe BKR sfafisfic as a measure of fluclualion of fhe error term remains frue 
in general. We have fhe following identify for compacfly supported lesl function^ 


S,{f) = -Jfix)s(^(0t + ^)dx + O{\/X{t)). (16) 

Clearly, fhis equation recovers (fT4b in fhe case of f = Xu since f is simply fhe difference of fhe della 
functions al fhe endpoinls. 

We now proceed to slate our resulls. Lef 0 < « < 1 and Xu\x)^ i = 1,2, denote fhe indicator 
functions of fhe intervals [0, u] and [—f/e,u] for some fixed e > 0, respectively. Lef 0(x) be a 
smoofh bump function supported on (—1/2, 1/2), and denofe 


K = — j <^{x)<^{y)\og\x — y\dxdy. (17) 

^ This equation does not appear in GS but follows easily from intermediate steps in their derivation of Proposition 3. 
In all of our applications, cf. 1211 and 1361 below, the constant in 0(l/A(f)) is of the order 0(| |/|| i) so that the 0(l/A(f)) 
term is negligible for our purposes due to 11/| | i/A (?) <C 1. A proof of is given in the appendix. 
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Define the £-rescaled bump function by 0e(jc) = 1 /e^{x/e), and let felix) be the smoothed indicator 
functions of the intervals [0, u] and [—!/£,«] given by the convolution of Xu\x) with 0e(^) 

/ijw - = “/ Xu\x-y)^{y/£)dy. (18) 

Clearly, (x) = 1 for ;k: G [£/2, n — £/2], (x) = 1 for x G [—1/£ + £/2, m — £/2], /g'j (x) = 0 for 

X > M + £/2 and x < —ejl, and /g^j(x) = 0 for x > m + £/2 and x < —lle — ejl. Moreover, 

/el\x) =<l>e{x)-(j)e{x-u), (19) 

fel\x) =^e{x+l/e)-(l)e{x-u). (20) 


Theorem 12. 1 1 applies to feji{x) for all m > £ > 0. Fix £ > 0 and denote the BKR statistic based on 
some constant times f^l (x) by sf ^ (/!,«,£) 


sl’\n,u,e) = 71 *</>£)( 7 ( 0 ) - 




log? 

27iX{t) 


(0 


*(/)g)(x)dx 


( 21 ) 


The meaning of the first statistic (t = 1) is that it counts zeroes in the interval [t,2l] over three asymp¬ 
totic scales. Specifically, over the interval of length u/X{t) — e/X{t), 


u £ 


<(Ot < 


£ 

l^0’ 


( 22 ) 


the zero is counted with the weight 1, whereas it is counted with a diminishing weight that is deter¬ 
mined by £ and the bump function over the boundary intervals of length e/X{t), 


E E 

y- w^rr\ < y^ 


2X{t) 2A(0’ 

£ U 

<(Ot < Y-^r^ + 


' X{t) 2A(0 A(0 2A(0' 

The second statistic (/ = 2) has the same interpretation except that instead of 


7- 

7- 


■-f ■ 


A(0 2A(0 

£ 


<cot < y- 


+ 


1 


2A(0 eA(0’ 


■-f ■ 


2X{t) EX{t) 


1 £ 
<(ot < y+ 


+ 


2A(0 eA(0’ 


(23) 

(24) 

(l22l) and (|2^ we have 

(25) 

(26) 


respectively, so that the middle interval has length u/X{t) — e/X{t) + l/eX{t), whereas the boundary 
intervals have length e/X{t) as before. Hence, the three scales are e/X{t), u/X{t), t, and they satisfy 
the asymptotic condition 


1 £ 
log? X{t) 


< 


u 

W) 


<Cl <?. 


(27) 


This means that t defines the global scale and ujXit) S> average spacing so that u/X{t) and e/X{t) 
are on the mesoscopic scale. 


*The choice of the constant 7ty/2fl, 0 < ^ <2 will be explained in Section 4, cf. Corollarv l4.2l 
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The meaning of the statistics sl’\lJ.,u,e) can be elucidated further by means of Theorem 12. 1 1 and 
(fT^ . Given the expressions for the derivatives in ([T^ and (l20l) . which show that the derivatives are 
smooth approximations of the difference of the delta functions at the endpoints, we can write 


St^\n,U,£) = j - u) - (i)E{x)) s{G)t + dx+ 0[\/ X{t)) , 


St^\n,U,£) = J [(l)e{x — u) — 0e(x+ l/s)) s(^(Ot + dx + 0[\/X{t)) 


SO that the statistics are smooth approximations of 




A(0 


u 


X{t) 


-S{(ot) 

-sicot- 


1 


£X{t)J. 


(28) 

(29) 

(30) 

(31) 


and can be interpreted as fluctuations of the smoothed error term over [cut, (ot + u/X{t)] and [cot — 
l/eA(t), (Ot + u/X{t)], respectively. Moreover, by Theorem 12. 1 1 and calculations in Corollaries 14.21 
and 14.41 the processes u —)• S^\ix,u,£), u E (0,1), converge in law in the limit t —oo to centered 
gaussian fields having the asymptotic covariances, for / = 1, 


! — ullage — tc + loglw — v| — log| m| — log|v|) + 0(e), if |m — v| S> £, 
—2/r(loge — fc —logluQ +0{e), if « = v, 

and, for i = 2, 

{—/i(31og£ — tc + loglu — v|) + 0(e), if |m — v| ^ e, 

— /i(41og£ — 2fc) +0(e), if M = V, 


(32) 


(33) 


thereby exhibiting logarithmic correlations. These covariances will be identified in Section 4 as corre¬ 
sponding to the centered gaussian free field and the gaussian free field plus an independent gaussian 
random variable, respectively. We see from (l32l) and (|33] ) that the asymptotic covariances of both 
statistics diverge as log£ in the limit £ —)■ 0. We are primarily interested in the statistical structure of 
the log |m — v| terms, which is hiding behind these divergences. We will remove them so as to reveal 
the underlying structure by introducing rescaling factors into the Mellin transforms as shown below. 

It should also be noted that the need for smoothing is necessitated by the singularity of our statis¬ 
tics in the limit £ —)• 0. Indeed, sf'\lJL,u,e) is not even defined in this limit, and Var[5'f^^(/i,M,0)] oc 
log(logf/A(t)) in Fujii’s theorem, cf. (fTSl) . is asymptotically different from Var[limSf^^(/x,M,£)] oc 

_ 7—^•oo 

— log(£) in (l32l) . This is explained by the difference between the formulas for the variance in (|4ll and 
(ITI), in fact, as shown in Lemma l4.101 the difference between the two for smoothed indicator functions 
is of the order 0(X(t)/elogt) and so becomes significant if one takes the £ —^ 0 limit first. To bypass 
this singularity, in what follows we will always take the t limit first and, when the two limits are taken 
simultaneously, e(t) will vary “slowly” enough to achieve he same end. 

Motivated by the proof of Theorem l2.2l in |[T5l . we are also interested in the same type of statistics 
that are based on f-dependent e(t). Let 


logt 


logt 


<C e(t) <C 1, for / = 1, 


1 


< e(t) <C 1 and for * = 2- 


(34) 

(35) 
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In other words, we replace e with e{l) in (l2TI) . resulting in the statistic^ 


logt 




^<Pe{t)){x)dx 


(36) 


These statistics also count zeroes in the interval [l,2t] over the three scales e{t)/X{t), u/X{t), t, as 
before, except that u/X{t) S> e{t)/X{t) average spacing, 


1 eit) u 

— <c —^ <c — 

logt X{t) ^ X{t) 


< 1 <t. 


(37) 


Remark We note that the asymptotic conditions in (l34l) and (l35l) have a simple interpretation of being 
“slow” decay conditions that require that the lengths of intervals, over which the zeroes are counted, 
go to zero, whereas the expected numbers of the zeroes counted go to infinity. The lengths of the mid¬ 
dle and boundary intervals of the first statistic satisfy 0(l/A(t)) and 0{e{t) , so that the corre¬ 
sponding expected numbers of zeroes are C?(logt/A(t)) and 0(^e{t) logt/A(t)) by (fT^ . respectively. 
Thus, the lengths go to zero by (l2]l and the expected numbers go to infinity by (l2]l and (l34l) . Similarly, 
the lengths of the middle and boundary intervals of the second statistic are 0{\/X{t) + l/e(t) A(t)) 
and C?(e(t)/A(t)) so that the above argument remains valid provided e{t)X{t) S> 1, hence (1351) . We 

finally nofe fhaf fhe inferprefafion of (ft,«) as flucfuafions of smoofhed error ferms and calculafions 
of limiting covariances in (l32l) and (l33l) remain valid provided (l34l ) and (1^ are salisfied as shown 
in Lemmas 14.91 and |4.101 i.e. fhe limiting fields are gaussian and have £(t)-dependenl asymptotic 
covariance given in (l3^ and (1331) wifh e = e{t). 


The sfafemenfs of our resulfs below require some familiarily wifh fhe double gamma funclion 
of Barnes (or fhe Alexeiewsky-Barnes G—function). We will give a brief summary of ifs definifion 
and properfies here and refer fhe reader to ftil, f7l for the original construction, to lf54]| for a modem 
treatment, and to BSl and |[50l for detailed reviews. One starts with the double zeta function 


1^2{s,w\t)= Y, iw + ki+k2r) " (38) 

ki,k2=0 

that is defined for Re ( 5 ) >2, Re(>v) > 0, and T > 0. If can be analyfically exfended to 5 E C wifh 
simple poles af 5 = 1 and s = 2. The double gamma funclion is Ihen defined as fhe exponential of fhe 
^-derivative of (^2(‘5 ,w|t) al 5 = 0, 

r2(w|T) = exp(^5^|^^Q(^2(‘y,w|T)^. (39) 

The resulling funclion can be analyfically exfended lo w E C having no zeroes and poles al 


w = — (ki-|-k2T), ^1,^2 E N. (40) 

Barnes gave an infinile producl formula for r 2 (>v | t), which in our normalization lakes on fhe form 

“ / 

r2'(w|T) = /("'i")w n '(i + - 


ki,k2=0 


W 


k\ -\- k2'V 


exp 


w 


■ + 


w 


ki k2T 2{k\-\-k2T') 


2 / ’ 


(41) 


®We note that Xii^\x) = X[-i/e{t) i<](x) becomes r-dependent so that ||;i;P^||i = 0(l/e(t)) and the 0(l/A(f)) term in 
GHi is negligible due to e{t)X(t) ;§> f in 














where P(w | t) is a polynomial in w of degree 2, and the prime indicates that the product is over all 
indices except k\=k 2 = 0. The double gamma function satisfies the functional equations 


^2{z\r) 
r2(z+l|T) 
r2(z|T) 
r2(z + T| t) 


VW’ 




r(^), 


(42) 

(43) 


where r(z) denotes Euler’s gamma function. An explicit integral representation of logr 2 (w| t) and 
additional infinite product representations of r 2 (w | t) can be found in Il48l . 

From now on, it is always assumed that A(t) and e{t) satisfy (jj)) and (l34l) or (1^ . respectively, 
and K is as in (fTTI) . Given 0 < p < 2, henceforth let 

T=->1. (44) 

P 

The results given below are stated for the Sp^(p,M,e) and Sf^\lJ.,u) statistics to avoid redundancy, 

(2\ 

for the same formulas apply to S, {ll,u,e) provided one simultaneously replaces 


^^ ^t^Xogeq^ 


(45) 

(46) 


and, similarly to (/!,«), provided one replaces <—>• and uses (1^ instead 

of (1^ . For clarity, this translation is shown explicitly in Conjectures 12.31 and 12.41 Given these 
preliminaries, our results are as follows. 


Conjecture 2.3 (Mellin-type transforms: weak version) Let Re{q) < i, then 


<?(<?+!) 


lim^M(loga-v)ii^ 

e->0 


limE 

/—>co 


= lim e‘ 

£-->0 


log f,^^hoge-K) 


•1 


limE 

t^oo 


gSl \n,u,e)^^ 


0 


= T?(27r)^r-‘?(l-l/T) 


r^(l -g + T| t) r2(-<? + T| t) r2(2-^ + 2T| t) 
r^(l + T|T) r2(T|T) r2(2-2(7 + 2 T|T) 


Let — (t+ 1)/2 < Re(^) < T, then 


(47) 

(48) 

(49) 


g(g+i) 


lim 

£-->0 


limE 

f—>co 


■?(■?+!) 


= lim 

£->0 


= T?(27r)^r-^(l-l/T) 


>1 


limE 

/—>oo 


r2(i+^+ t|t) r2(i — qx\x) r2(—<7+'f| x) r2(2+^+2T'| x) 


0 


r2(l + 2<7 + T I t) r2(l + T|T) r2(T|T) r2(2 + 2T'|T') 


(50) 

(51) 

(52) 
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Conjecture 2.4 (Mellin-type transforms: strong version) Let Re(^) < t, then 


limgM{iogeW-'c)^E 

/—yoo 


-M. 




lim e^ ® e^ (log e(0 -£ 

f—>co 




e^Q 




limE 

_f—>oo 


-M, 


,S?\n,u,e)du 


= lim '°g ('"g ® ^ * 

e->0 



limE 



.f —>00 

- \Jo '' -1 - 


(53) 


(54) 


Lef —(t+ 1)/2 < Re(^) < T, f/ien 




= lime^(^°g'' '"('A’ 

limE 


f —>00 

L \Jo ' 

e->0 

Y —>00 

L Vjo ' J J 


(55) 


lim e ^® hWe^ (log e(0 - 'f) £ 

f—>oo 






= lim 

e->0 


ic 


lim E 

f—>oo 




(56) 


The following corollaries can be formulated with either the double or single limits as in Conjectures 
[23] and [231 We will state them with the single limit for notational simplicity. All corollaries apply to 
both statistics, cf. ([45l) and ([46l) for the translation. Note that in Corollaries 12.51 [231 and l2.7l we have 
q = n, —n, N, respectively, so that the multiplier in (1451) and the scaling factor in ([4^ are adjusted 
accordingly, when translated for the second statistic. 


Corollary 2.5 (Positive integral moments and Selberg integral) Let n E N such that n < t. 
lime^*''°®®*'^^^'^^ ^2 *E (^J = J l^,-— (57) 


KJ 


[ 0 , 1 ]" 

«_l r(l-(;^+l)/T) r^(i-VT) 

M r(i-iA) r( 2 -(n+^-i)/T)’^ ^ 




n n 


-- / ^dsi---dsn, 

[ 0 , 1 ]" 


(59) 


r(i - (it+ i)/t) r(i + 2 n/T-it/'!^)r(i - 


n— 1 

n r(i-iA) r(2+n/T-(^-i)/T) 

(60) 


Corollary 2.6 (Negative integral moments) Let n E N. 


lim et‘('°g^(0-'^)A-^£ 

f —>00 




r(2+ (n + 2 + A:)/T)r(l - 1 /t) 
LI r^\ + ik+\)/r)r{l+k/r) 


(61) 
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Let n gN such that n <[t + l)/2. 


/—>00 




n— 1 


=n 


T{l+{-n + l + k)/z)T{\-\/z) 

T{l-2n/T+{k+l)/T)T{\ + {k+\)/x)T{l+k/x) ■ 

(62) 


Corollary 2.7 (Joint integral moments) Let m G N an<i denote N = n + m,N <2 /t. Let l\ and h 
be non-overlapping subintervals of the unit interval. Then, 


/—>00 




/—yoo 





A? 


— ^dsi---dsN, 
i<j 

(63) 





mn 



^dsi ■ ■ -ds;^. 
(64) 


Corollary 2.8 (Asymptotic expansions) Given q gC, then in the limit /i —)• 0 we have the asymp¬ 
totic expansion^^ 


limeM(iogeW-'^)^E 

t —^00 





00 



/xy+i 1 

2 ) 7TT 


C(r+l)x 


Br+2 (^ + 1) + 2Br+2 (^) ~ 35^+2 
r + 2 


-9] + (C(<-+l)-l)[ 


Br+2(g- 1) - Bf+2(2^ - 1] 


r + 2 


(65) 


lim e 

/—>oo 




~E 




■exp(£ 


r=0 


^uy+i 1 


r + 1 - 


^(r I I 


r + 2 


®r+2( <?) ®r+2( 2^) 26^+2 -| ^ ^'jf®r+2( <J 1) ^r+li, 1)- 

+--,] + (C(r+i)-i)[- — - 

Corollary 2.9 (Covariance structure) Let 0 < < 52 < 1- Then, in the limit A — )• 0, 


( 66 ) 


ii+A 


i2+A 


lim 

/—yoo 


/x(log£(f) — jc)+ Cov^log j du, log j 




S2 


il+A 


i2+A 


lim 

/—>00 


/x(31oge(f) — jc)+ Cov^log y e‘^'log y e'^' 


S2 


= -/l(log| 5 l - 52 I -log| 5 i|- 
-log|^ 2 |)+ 0 (A), (67) 

= - plog\si - 52 I + C?(A). 

( 68 ) 


The next two results deal with the probabilistic structure that is underlying the Mellin-type transform 
in (|4^ . 

As usual, ^ (i) denotes the Riemann zeta function and B„ (s) the nth Bernoulli polynomial. Note that (1) never enters 
any of the final formulas as the coefficient it multiplies is identically zero. 
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Corollary 2.10 (Non-central limit) The limit in (|4^ is the Mellin transform of a positive probability 
distribution. Call it M^. Then, for Re{q) < T, 


>oo 




0 ) 






(69) 


logM^ is infinitely divisible on the real line having the Levy-Khinchine decomposition logE[M^] = 
qm{p) + jq^o^{p) + /r\{o}(^^“ ~ 1 —qu/{\ + u^))d.Jf^{u) for some xn(p) G M and the following 
gaussian component and spectral function 




4 log 2 


(u = 


/ 


,-41+ t )/2 


dx 


(e^-l)(e^^-1) (W2-l)(t4V2_i)J X 


(70) 

(71) 


for u>0, and (u) = 0/or u<0. is a product of a lognormal, Frechet and independent Barnes 

beta random variables!^ 


Corollary 2.11 (Multifractality) Let 0 < 5 < 1. Then, for Re(^) < T, the limit 


lim 

/—>oo 





= eW(^)] 


(72) 


is the Mellin transform of a probability distribution, call it ( 5 ), which satisfies the multifractal 

lavE 

M^{s)^'^= stxp{a,)M^, (73) 

where Lis ci gaussian random variable with the mean {p/2)\ogs and variance —plogs that is 
independent ofM^ and is as in Corollarv \2.10\ This law is understood as the equality of random 
variables in law at fixed s < 1. In particular, M^{s) is also log-infnitely divisible. 


Corollary 2.12 (Multiscaling) Let 0 < 5 < 1. Then, for 0 < Re(^) < t, 


=const{p,q)s‘>-^^‘^"-‘^\ 


c.(l) 


limgM{ioge(f)-r)2i^E ( =const{p,q)s‘i+^^‘>"+‘^\ 

LVJo 7 J 

i.e. these Mellin-type transforms are multiscaling as functions of s. 


«(')/• 


(74) 

(75) 


The rationale for considering two separate transforms in Conjectures 12.31 and 12.41 is that they are 
complementary. The transforms in (1471) and (l5^ are not, strictly speaking, rescaled Mellin transforms 
because of the multiplier. This multiplier is introduced to obtain the Mellin transform of a 
probability distribution on the right-hand side of these equations. If we drop this multiplier, we loose 
this probabilistic interpretation but gain a bona fide rescaled Mellin transforms of fg e^' 

"This decomposition is shown in detail in Section 3, cf. Theorem l3.4l 

'^Also referred to in the literature as stochastic self-similarity or scale-consistent continuous dilation invariance. It 
determines the law of (/,■?< 1, from the law of in l l69t . 
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and /o' on the left-hand side so that the limit fits into the general theory of mod-gaussian 

convergence, cf. ||33l. The interest in introducing the second statistic sf'\lJL,u,e) and sf'\lJL,u) is 
precisely to eliminate the need to introduce the multiplier so that the transforms in (|4^ and 
(l54l) are both bona fide rescaled Mellin transforms and give the Mellin transform of a probability 
distribution on the right-hand side of these equations. We note that the right-hand sides of both (|4^ 
and (l52l ) are special cases of the Mellin transform of the Selberg integral distribution that is given in 
(II171) and (I135I) . The meaning of the conditions Re(^) < T and —(t-|- 1)/2 < Re(^) < T is that the 
right-hand sides of (|4^ and (l52l ) are analytic and zero-free over these regions. As will be explained 
in Section 4, one can insert the (1 — u)^- prefactor into the exponential functionals on the left-hand 
sides of Wh - (|5^ so as to obtain the general Mellin transform of the Selberg integral distribution on 
the right-hand side, cf. Theorem 14. 5 1 for details. We chose not to state the most general case here to 
avoid unnecessary complexity. 

Remark A random matrix analogue of the Bourgade-Kuan-Rodgers theorems was established in 
1271. The corresponding linear statistic is an appropriately centered log-absolute value of the charac¬ 
teristic polynomial of a suitably scaled GUE matrix. The eigenvalues of the matrix replace the roots 
and its dimension N replaces logt. It is shown in 1^ that the A —)• oo limit of the GUE statistic is 
gaussian, and the limiting covariance can be showrl^to be equivalent to the Bourgade-Kuan-Rodgers 
covariance. As it will become clear in Section 4, our results apply to any linear statistic that is 
asymptotically gaussian with the Bourgade-Kuan-Rodgers covariance. Thus, under the correspon¬ 
dence N ~ logt, they apply equally to the GUE statistic. The significance of this observation is that 
the convergence of our statistics to the gaussian free field or its centered version provides a theoretical 
explanation for why our conjecture about the zeroes can be approached directly from the GUE side. 

Remark Our conjecture has only GUE terms and no arithmetic corrections, contrary to most of the 
conjectures that we mentioned in the Introduction. The reason for this is the choice of the mesoscopic 
scale 1 <C A(t). As Eujii |[24l shows in his analysis of Berry’s conjecture, which in particular deals 
with the second moment of the fluctuation of the error term, cf. (fT4b . the GUE term in Berry’s formula 
dominates arithmetic corrections precisely under the condition 1 <C A(f). As our conjecture involves 
essentially the same statistic, except for smoothing, we expect that in our case this condition achieves 
the same effect of dominating arithmetic terms. 

Remark The interest in the strong version of our conjecture is that it contains information about the 
statistical distribution of the zeroes at large but finite t, whereas the weak version only describes the 
distribution at t = oo. Indeed, as it will be explained in Section 4, the strong version is equivalent to the 
weak version provided the statistics s\'^ {}X,u) converge to their gaussian limits faster than 1/| log£(t)|, 
in the sense of the rate of convergence of the variance to its asymptotic value, and this is expected 
to be the case due to the asymptotic condition in (l34l) . Moreover, as the strong conjecture fits into 
the general framework of mod-gaussian convergence, the results of lITTI and |[39l and the explicit 
knowledge of our limiting functions make it possible to quantify the normality zone, i.e. the scale up 
to which the tails of our exponential functionals are normal, and the breaking of symmetry near the 
edges of the normality zone thereby quantifying precise deviations at large t. The actual computation 
of these quantities is left for future research. 

'^Nickolas Simm, private communication. 
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3 A Review of the Limit Lognormal Measure and Selberg Integral Dis¬ 
tribution 

In this section we will give a self-contained review of the limit lognormal measure of Mandelbrot- 
Bacry-Muzy on the unit interval and of the Selberg integral probability distribution mainly following 
our earlier presentations in P3I . 1481 . and l50l . We will indicate what is known and what is conjec¬ 
tured and refer the reader to appropriate original publications for the proofs. 

The limit lognormal measure (also known as lognormal multiplicative chaos) is defined as the 
exponential functional of the gaussian free field. Lef co^x,e{^^) be a sfafionary gaussian process in 5 , 
whose mean and covariance are functions of fhree paramefers /i > 0, L > 0, and e > 0. We consider 


fhe random measure 

b 

= (76) 

a 

The mean and covariance of defined fo be, cf. flTl . 

E[«^,L,r(0] =-f (l+log^), (77) 

Cov [a)^x,e(f), =Atlog <L, (78) 

Cov [a)^,L,e(0: £(■«)] = M n-f log^ j , |f-5| < £, (79) 


and covariance is zero in fhe remaining case of |t — i'| > L. Thus, £ is used as a fruncafion scale. L 
is fhe decorrelafion lengfh of fhe process, /r is fhe infermiffency paramefer (also known as inverse 
femperafure in fhe physics liferafure). The fwo key properties of fhis consfrucfion are, firsf, fhaf 

E[a)^,L,e(t)] =-^Var[a)^,L,e(t)] (80) 

so fhaf 

E[exp(a)^,z._e(5))] = 1 (81) 

and, second, fhaf 

Var[a)^,L,e(t)] oc-loge (82) 

is logarifhmically divergenf as £ —;■ 0. The firsf properfy is essenfial for convergence, fhe second is 
responsible for mulfifracfalify, and bofh are originally due fo Mandelbrof l37l . 

The inferesf in fhe limif lognormal consfrucfion sfems from fhe £ —>• 0 limif. If is clear fhaf fhe 
£ —)• 0 limif of does nof exisf as a stochastic process (fhis limiting “process” is known as fhe 

gaussian free field). Remarkably, using fhe fheory of T-marfingales developed by Kahane 1341 . fhe 
work of Barral and Mandelbrof f91 on log-Poisson cascades, and conical consfrucfions of Rajpuf and 
Rosinski lISTI and Marsan and Schmiff |[55l . Bacry and Muzy O showed fhaf M^_i,e(<if) converges 
weakly (as a measure on M+) a.s. to a non-frivial random limit measure 

= ( 83 ) 
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provided 0 < At < 2 , and the limit is stationary M^ L{t,t + t) '"=* t) and non-degenerate in 

the sense that E[M^ ^(0, t)] =t. We will denote the total random mass of [0,L] by 

= (84) 


It is shown in O that for ^ > 0 we have 

2 


q< 


E [M^ l] < ^ ~ = oo. 


(85) 


The fundamental property of the limit measure is that it is multifractal. This can be understood 
at several levels, for our purposes, this means that its total mass exhibits stochastic self-similarity, 
also known as continuous dilation invariance, as first established in Q. Given y < 1, let Q.y denote a 
gaussian random variable that is independent of the process (O^y^ei^) such that 


= ^jUlogr, (86) 

Var[fly] = — Atlogy. (87) 

Then, there hold the following invariance^ 

COii,L,reir^), ( 88 ) 


that are understood as equalities in law of stochastic processes in s on the interval [0,L] at fixed £, L, 
and 0 < 7 < 1. In (l90l ) the superscripts denote independent copies of the free field, e sfands for fhe 
base of fhe natural logarithm, and ^ < 7 < 1. The truncation scale invariance in (1^ implies the 
multifractal law of the limit measurJ^ for f < L 


1-r //^ \ inlaw ^ 0.1 

4T^,l(0, 0 = -e 


(91) 


which implies that the moments of the total mass obey for 0 < <7 < 2 /At the multiscaling law 


E[M^x(0,fr] 




(92) 


as a function of t < L. The decorrelation length invariance in (IMl) implies that the dependence of the 
total mass on L is trivial!^ 


, , inlaw rut 

— LM^i, 


(93) 


*^The first invariance was discovered in and later generalized in ED- We discovered the other two invariances in |43I 
and developed a general theory of such invariances in (46|. 

*^It must be emphasized that this equality is strictly in law, that is, Sly is not a stochastic process, i.e. Sly and for 
y / / are not defined on the same probability space. In particular, l l91t determines the distribution of l) in terms 

of the law of the total mass but says nothing about the latter or their joint distribution. 

*®This invariance also determines how the law of the total mass behaves under a particular change of probability measure, 
cf. (46) for details. 
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so that we can restrict ourselves to L = 1 without a loss of generality. Henceforth, 


L = l, 


( 94 ) 


and L is dropped from all subsequent formulas. Finally, the significance of the intermittency invari¬ 
ance in (l90l ) is that it gives the rule of intermittency differentiation and effectively determines the law 
of the total mass, cf. Theorem 13. II and the discussion following it below. 

The law of the total mass can be reformulated as a non-central limit problem. Let us break up the 
unit interval into the subintervals of length e so that Sj = je, (Oj = and = 1. It is shown 

in lO that the total mass can be approximated as 

N-l 

e ^ M^as e ^ 0. (95) 


The essence of this result is that the limit is not affected by one’s truncation of covariance so long 
as dSnil holds. This representation is quite useful in calculations. In particular, it is easy to see 
that it implies an important relationship between the moments of the generalized total maso of the 
limit measure and a class of (generalized) Selberg integrals, originally due to lO. Let (p{s) be an 
appropriately chosen test function and I a subinterval of the unit interval, then 


E 


(p{s)M^{ds)y 



(96) 


More generally, the same type of result holds for any finite number of subintervals of the unit interval. 
For example, the joint {n,m) moment is given by a generalized Selberg integral of dimension n-\-m 


E 


[J <?’i(5)M^(d5)) 

h 


[J 9i{s)Myds)^ 

h 



n n+m n-\-m 


n '(’2(5,on 

i—\ i—n+l k<p 


tp 


4 1 ^dS\---dSn+m- 
(97) 


It can be shown as a corollary of this equation, cf. f3i], that the covariance structure of the total mass 
is logarithmic. Given 0 < ^ < 1, 


i+Ai Ai 

Cov^log J M^{du), log J M^ldu^j =—iJ.logs + 0{As), (98) 

i 0 

which indicates that the mass of non-overlapping subintervals of the unit interval exhibits strong 
stochastic dependence. In the special case of q){s) = a^'( 1 we have an explicit formula for 

moments of order n < 2/jj., as was first pointed out in ||4l, that is given by the classical Selberg 
integral, cf. Chapter 4 of 1221 for a modem treatment and |[5^ for the original derivation. 


E 



-sf^Myds)y 


r(l - (k + l)At/ 2 )r(l + Ai - kn/ 2 )r{l +X 2 - k^/ 2 ) 

r(l -At/2)r(2 + Xi+X 2 -in + k- l)At/2) 


a slight abuse of terminology, we refer to any integral of the form Jg* (p{s)Mji (ds) as the generalized total mass. 
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The law of /J 5^* (1 — is not known rigorously, even for Ai = A 2 = 0, i.e. the total mass. 

It is possible to derive it heuristically so as to formulate a precise conjecture about it as follows. 
Consider the expectation of a general functional of the limit lognormal measure 




F 





( 100 ) 


The integration with respect to M^{ds) is understood in the sense of e —)• 0 limit so that v(/i, (p, f, F) = 


limve(At, q>,f,F) and v^in, q),f,F) = E F\^Jq 
Also, let g(5i, S 2 ) be defined by 

■^ 2 ) = -logl^'l -52|- 


with g{ds) as in (1761 ). 

( 101 ) 


Finally, we will use [0, 1]^ to denote the k—dimensional unit interval [0, 1] x • • • x [0, 1]. Then, we 
have the following rule of intermittency differentiation, cf. Il4^ and Il48l for derivations and iHTI for 
an extension to the joint distribution of the mass of multiple subintervals of the unit interval. 


Theorem 3.1 (Intermittency differentiation) The expectation v(/r, (p, /, F) is invariant under in¬ 
termittency differentiation and satisfies 


^v{p,(p,f,F)= I v{p,(p,f + gi;s),F^^^)e^f^-^^fis)(p{s)ds+ 

[ 0 , 1 ] 

/ ^{l^^V^f + si-,^i)+g{-,^2),F^^^)e^i^^'''>^^^'^'^'^^^'^’'^'^)g{si,S2)q>{si)q){s2)dsids2. ( 102 ) 
[o,i]2 

The mathematical content of (11021) is that differentiation with respect to the intermittency parameter p 
is equivalent to a combination of two functional shifts induced by the g function. This differentiation 
rule is nonlocal as it involves the entire path of the process s —)• M^{0, 5 ), 5 € (0, 1). It is clear that 
both terms in (11021) are of the same functional form as the original functional in (1 1001) so that Theorem 
13. H allows us to compute derivatives of all orders. There results the following formal expansion with 
some coefficients that are independent of F. Let 


„ / i 

Si{p,(p)= I ]^(p(5;) 

(=1 i<j 


[ 0 , 1 ]' 


A 

X = 


1 

j (p{s)ds. 


(103) 


(104) 


0 

Then, we obtain the formal intermittency expansion 

/o 


E 


00 rt r 2n 


F[l tp{s)Mffds)) =F{x)+'£^ £f«(x)//„,,(.p) 


n=\ "k=2 


(105) 


The expansion coefficients Hn^k{(p) are given by the binomial transform of the derivatives of the 
positive integral moments 


Hn,k{tp') — 




^"5/ , 

dff 


[/j=0) 


(106) 
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and satisfy the identity 


Hn.k{(p) = 0 Vfc > 2n. 


(107) 


Remark The last equation and Theorem 13.11 say that the intermittency expansion in (11051 ) is an 
exactly renormalized expansion in the centered moments of (p{s)dM^{s). Indeed, we have the 
identity 





(p{s)M^{ds) 



= k\Hn,k{(p)- 


(108) 


From now on we will focus on 




(109) 


which corresponds to the full Selberg integral. Clearly, we have 

_ r(i+Ai)r(i+A2) 
r(2+Ai + A2) 


( 110 ) 


The moments S/(/i, q)) = ^/(/i, Ai, A 2 ) are given by Selberg’s product formula in (l99l) . By expand¬ 
ing log5'/(/i, Ai, A 2 ) in powers of jj. near zero and computing the resulting Hn^k{(p) coefficients, we 
derived in |(45]1 and 1(48 1 the following expansion for the Mellin transform in terms of Bernoulli poly¬ 
nomials and values of the Hurwitz zeta function at the integers 


E 


f>r(9, Al, A 2 ) = Tt [(«('■+ 1,1 + A,) + C(|-+1, 1 +A 2 )) 

/ Br+2(^+ 1) —Br+2\ 


( 111 ) 


-C(^+l)9 + C(^+l) 

— ^(r 1,2 + X 2 ) 


r + 2 J 

Br+2{2q-l)-Br+2iq-i) 
r + 2 


( 112 ) 


The series in (11111) is generally divergenf^ and interpreted as the asymptotic expansion of the Mellin 
transform in the limit /r —)• 0. The intermittency expansion of the Mellin transform implies a similar 
expansion of the general transform of logM^(Ai, A 2 ), where we introduced the normalized distribu¬ 
tion 

r(2 + Ai+A2) fi 


M^(Ai, A2) — 




r(i -|- Ai)r(i -|- A2) Jo 

Then, given constants a and s and a smooth function F{s), the intermittency expansion is 


E 


n—0 


f ( 5 -halogM^(Ai, A2)) = Ai, A2)^, 


Foia, s, Ai, A2) =F{s), 


(113) 

(114) 

(115) 


Fn+i{a, s, Ai, A2) = Yj Ai, A2^F;-r(a:, 5 , Ai, A2). ( 116 ) 
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It is convergent for ranges of integral q as shown in 1451 . 
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The expansion is thus obtained by replacing q with ad/ds in the solution for the Mellin transform in 
(11111) . We note parenthetically that the general transform is particularly interesting in the special case 
of apurely imaginary a, in which case the operator F{s) —)■ E F(^s A 2 )) is formally 
self-adjoint. 

The calculation of the intermittency expansion of the Mellin transform naturally poses the prob¬ 
lem of constructing a positive probability distribution such that its positive integral moments are given 
by Selberg’s formula, cf. (|9^ . and the asymptotic expansion of its Mellin transform coincides with 
the intermittency expansion in (II111) . Equivalently, one wants to construct a meromorphic function 
I /I, Ai, A 2 ) that (1) recovers Selberg’s formula for positive integral ^ < 2//r, (2) is the Mellin 
transform of a probability distribution for Re(^) < 2//i as long as 0 < /X <2, and (3) has the asymp¬ 
totic expansion in the limit /r —)■ 0 that is given in (11111) . Such a function can be naturally thought of 
as an analytic continuation of the Selberg integral as a function of its dimension to the complex plane. 

We will describe an analytic and a probabilistic solution to this problem that we first found in | 
in the special case of Ai = A 2 = 0 and then in 1481 in general!^ Define 


m{q\n,Xuh) = T'(27rrr-'?(1 - 1 /t) 


r 2 (l —<7-|-t(1 -|-Ai) I t) r 2 (i — q-\- -|- A 2 ) I t) 


r2(i -|- t(i -i- Ai) I t) 
r2(~^-|-T|T) r2(2 — ^t( 2-|-Ai-|-A 2 ) I t) 
r2('f|T) r2(2 — 2 ^-|-t( 2-|-Ai-|-A 2 ) I t) 


r2(i -i- t(i -i-A 2 ) I t) 


(117) 


Theorem 3.2 (Existence of the Selherg integral distribution) Let 0 < /x < 2, A, > —lijl-, and z = 
2//I. Then, Tl{q | /I, Ai, A 2 ) is the Mellin transform of a probability distribution on {0,oo)for'Re.{q) < 
T, and its moments satisfy 


min IP,mm = n ni-(^+i)A)r(i+A,-V.)r(i+A 2 ^/T) ^ ^ ^ ^ ^ 

to r(i-i/T)r(2+Ai+A2-(n+^-i)/T) ’ - 

. , , , ^ r (2 + Ai+A 2 + (n + 2 + )t)/T)r(l-l/T) 


r(i -|-Ai -|- ik-\- i)/T)r(i -1-A2-1- (^-i- i)/T)r(i 

The function iXftiq \ p, Xi, X 2 ) satisfies the functional equations 


miq\p,Xi,X2) = miq-i\p,Xi,X2) 


r(i - ^/T)r(2+Ai+A2 - - 2 )/t) 


r(i-i/T) 

r(i-i-Ai — (^—i)/T)r(i-i-A2 — (^—i)/t) 

r(2 + Ai + A 2 - {2q - 2)/T)r(2 + Ai + A 2 - (Iq - 3)/t) ' 


(118) 

(119) 


( 120 ) 


m{q\p,XuM)=miq-T\p,Xi,?i2)ri2ny-M^(^\-^'^r{T-q)ir 

r((i+Ai)T-(<?-i))r((i+A2)T-(^-i)) r((2+Ai+A2)T-(^-2)) 
r ((2 + Ai+A2)T-(2<?-2)) r ((3 + Ai+A2)T-(2<7-2))’ 

( 121 ) 

*®The general case was first considered by Fyodorov et. at. 1281 . who gave an equivalent expression for the right-hand 
side of dl 17b and so recovered the positive integral moments without proving analytically that their formula corresponds 
to the Mellin transform of a probability distribution or matching the asymptotic expansion, i.e. solved (1) only. Instead, 
they used Selberg’s formula to deduce the functional equation in ( 1120b for positive integral q, conjectured that it holds for 
complex q, and then found a meromorphic function that satisfies ( inob . 
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The function log5[R(^ | /x, Ai, A 2 ) has the asymptotic expansion as p ^ +0 


logm{q\p,XuX 2 ) ~^log 


r(i+Ai)r(i+A2) 


r(2+Aj + A2) 

+ {l^{r+i-,i+Xi) + i;{r+l,l 
r + 2 / 




Br+2{2q - 1) - Br+2{q - f 
r + 2 


(122) 


The structure of the corresponding probability distribution, which we denote by most 

naturally explained using the theory of Barnes beta distributions that we developed in ||48l, For 
simplicity, we restrict ourselves here to the special case of the distribution of type ( 2 , 2 ) and refer the 
reader to BOll . ll50l for the general case. 


Theorem 3.3 (Barnes beta distribution of type (2,2)) Let bo,bi,b 2 >0 and T > 0. Define 

T2{b{) + b2\'c) T2{q + bQ + bi+b2\'c) 


B{q\T,b) = 


2^T2{q + bQ\z) T2(bQ + bi\z) 


r2(Z7o|T) T2{q + bQ + bY\'c)T2{q + bQ + b2\'c) T2{bo+b\+b2\'c) 

^ [ bo + n\+n2't: q + bo + bi+ni+n2't: q + bQ + b2 + ni+n2'C 


nin 2 =tfid + bQ + nx+n2T bo + bi+nY+n2T 
bo + bi +b2 + ni +n2T 


bo + b2 + ni+n2T 


(123) 


(124) 

q + bo + bi+b2 + ni +n2'ci 
Then, ri(q\T,b) is the Mellin transform of a probability distribution on {0,1). Denote it by l5{T,b). 

E [i3 (T, b)‘^] =r]{q\r,b), Re(q) > -bo- (125) 


The distribution —logl5{T,b) is absolutely continuous and infinitely divisible on (0,oo) and has the 
Levy-Khinchine decomposition 


E 


exp{qlogp{T:,b)) 


= exp 


{e-’“i-l)e-’’°’‘ 


{I - e-^'^){l - dx 
(l-r’^^)(l-e-«) ~ 


,R&{q)>-bo. (126) 


We can now describe the probabilistic structure of the Selberg integral distribution. Let T > 1 and 
define a lognormal random variable 


L = exp(./L(0, 41og2/T)), 


(127) 


where ^(0, 41og2/T) denotes a zero-mean gaussian with variance 41og2/T, and a Frechet variable 
Y having density Ty^'^^exp(— dy, y > 0, so that its Mellin transform is 

£[7^?] =r(l-|), Re(^)<T, (128) 

and logT is infinitely divisible. Given A,- > —1 /t, letXi,X 2 ,X 3 have the p^^{T,b) distribution with 
the parameter^ 

Xi ^j 82 : 2 '(g^o = 1 + t + tAi,(7i =t(A2-Ai)/2,(72 = t(A2-Ai)/2), (129) 

X 2 ^ = 1 + T + t(Ai + A 2 )/ 2 , bi = 1/2, b 2 = T/ 2 ), (130) 

X 3 = 1^2 2 = 1 + T, (>1 = (1 + T + tAi + 1 X 2 )(2, f ?2 = (1 + T+ tAi + 1 X 2 )I2)j. (131) 

Without loss of generality, A 2 > Ai. If A 2 = Ai, then Xi = 1, otherwise the parameters of Xi satisfy Theorem 1 3.3 1 
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Theorem 3.4 (Structure of the Selherg integral distribution) Let z = 2/jj.. decomposes 

into independent factors, 

- \/zY^LX^XxX^Y. (132) 


In particular, ^ogM(^^ Xi,k 2 ) absolutely continuous and infinitely divisible. Its Levy-Khinchine de¬ 
composition logim{q\p,Xi,X 2 )=qm{p)-\-Y^O^{p)-\-Y\^^Q-^ [e‘^'‘-l-qu/{l+u^))d.y£{^^^XuX 2 ){u) 
is 


o^{p) = 


4 log 2 


•^(M,Ai,A2)(“) = - J 


-I- _i_g-xTX2 _j_g-;c(l+T(l+Ai+A2))^ 
(e^-l)(e^T_i) 


(133) 


3:(1+t(1+Ai+A2))/2 - 
(f4/2- l)(f4V2- 1). 



for u > 0, and .^(^,A|,A 2 )(“) ~ ^ foz u <0, and some constant m{p) G M. The Stieltjes moment 
problems for both positive, cf (11181) . and negative, cf (II191) . moments of M(^^Xi,X 2 ) are indeterminate. 

We finally note that the Mellin transform of the Selherg integral distribution has a remarkable factor¬ 
ization, which extends Selberg’s finite product of gamma factors to an infinite product. 


Theorem 3.5 (Factorization of the Mellin transform) 


Tl{q\p,Xi,X2) 


^T'^r(l — ^/T)r(2 —2 ^-|-t( 1-|-Ai-I-A 2 )) 2^ r(l — 

“ R(l-l/T)r( 2 -^ + T(l+Ai+A2)) IJ, T{\+mz) 

r(l —q-\-zX\ -\-mz) r(l — ^TA2-|-mr) r (2 —^-|-t(Ai -I-A2) -\-mz) 
T(^1ZXimz^ r^l-f TA 2 “h r^2 — 2qz(^Xi^ 2 ^-\-mz^ 


(135) 


We conjectured in Il45l in the special case of Ai = A 2 = 0 and in Il48l in general that the Sel- 
berg integral distribution coincides with the law of the generalized total mass of the limit lognormal 
measure. 


Conjecture 3.6 Let p G (0,2) and A,- > —pl2, then 

j\^'{\-sYMYds) "=M(^^XuX2y (136) 

The rationale for this conjecture is that we constructed a family of probability distributions parame¬ 
terized by /I, Ai, and A 2 having the properties that (1) its moments match the moments of the (gener¬ 
alized) total mass Jq (1 — s)^^ (ds ), i.e. the Selherg integral and (2) the asymptotic expansion of 
its Mellin transform coincides with the intermittency expansion of the total mass. It is finally worth 
pointing out that the restriction A; > —p /2 is artificial and only imposed so that the Mellin transform 
can be defined over Re(^) <2jp. It is not difficult to see from (|9TI ) that 

+ (137) 

log! 2 

in the limit 1 —)• 0 so that the integral Y 5 ^' (1 — (ds) is actually defined for A,- -|- /r /2 > — 1 and 

its Mellin transform for Re(^ 7 ) < min{2//r, 1 -|- (1 +Xi)2/p}. 
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4 Calculations 


In this section we will give heuristic derivations of Conjectures 12.31 and 12.41 and their corollaries 
assuming Coniecture l3.6l to be true. The derivations are based on certain exact calculations, combined 
with a key assumption in the case of Conjecture 12.41 that is explained below. 

The main idea of the derivation of Conjecture I2.3l is that the 5^ (ft,£) and Sp^ (ft,£) statistics 
that we introduced in (|2TI) in Section 2 converge in law, up to 0{e) terms, to (modifications of) 
the processes ft)^ e(M) — (O^ g{0) and respectively, in the limit t ^ where (Op^^e{u) is the 

approximation of the gaussian free field that we defined in (iTTll -dT^ (recall L = 1). 

Lemma 4.1 (Main lemma for Let be as in (fT^ . Xu\x) = X[o,u]{^): 

/i,M W = {Xu^<^e){x) = ^J Xu\x-y)(j){y/e)dy, (138) 

K be as in dnii, and the scalar product be defined as in ©. Then, given 0 < u,v < I, in the limit 

£ —;■ 0 , 

(/PP/i.v) = -^(log£- K- + log|M-v| -login] -log|v|) +0(£), i/ln-vl >£, (139) 

= -^(loge-K--log|M|)+0(£), ifu = v. (140) 


Proof Recall that /PP = (x) — ^g(x — n), where ^g(x) = l/e^{x/e) is the rescaled bump function, 

cf. ([T^ . Then, 


/a t M I \n J log k +0 £p, if Z >£, 

/ (/)g(x)log|x-z|dx= (141) 

J I log£ + /0(x)log X- I dx, if |z| < £. 


It follows 


log jyj - logjn -y| + 0{e^) if |y|, ly- m| > £, 

log£ + /0(x)log |x- ||r/x-log|M-y| +0(£^) if jyj < £, 


J (0g(x)-0g(x-M))log |x-y|<ix = 

log jyj - log £ - / 0 (x) log |x - I ^i^x + 0{e^) 

if |y — m| < £. Now, by the definition of the scalar product, 

= -^ J ((^e(y)-0e(y-v))(0g(x)-(/)e(x-M))log|x-y|dxr/y, 


(142) 


= -^/(0eW-(/>g(x-M))log|x-y|r/. 

+ ^/ My) (^e(x)-^g(x-n))log|x-(y + v)|r/. 


dy 


dy. 


(143) 


The resulting integrals are all of the functional form that we treated in (11411) and (11421) so the result 
follows. I 


22 











Corollary 4.2 (Covariance structure of Let the statistic Si^\lJ.,u,e) be as in (l2T]) . 

A ^ /T: 7 rv 


S} \n,u,e) = 7rv^[L/eM7(0) - 2^x{t) 




(144) 


Then, in the limit f —>• oo the process u —>• Si^\p,u,e), u G (0,1), converges in law to the centered 
gaussian field S^^\p,u,s) having the asymptotic covariance 

C„v(sI')(M,»,E).5'‘)(f,.v,E))=| -Wloge-^ + logl«-v|-log|„|-logH),./|.-v|»e, 
V I -2At(log£-JC-log|M|),//■ m=v, 

+ 0{e). (145) 

Proof This follows from Theorem 12. II and the calculation of the scalar product in Lemma |4~T] | 


Lemma 4.3 (Main lemma for (/!,«,£)) Let f'fu[x) be as in (fT^ . Xu^x) = X[-i/e,u]{x)- Then, 
given 0 < M,v < 1, m the limit £ —)• 0, 


f(2)/ 


,(2)/ 


(/S,/S) = -^(31oge-K- + log|M-v|)+0(£), //■|m-v| >£, 
(/i2,/S) = -^(41oge-2K:)+0(£), ifu = v. 

Proof The calculation is the same as in the proof of Lemma |4~T] and will be omitted. | 


(146) 

(147) 


Corollary 4.4 (Covariance structure of sf'\p,u,e)) Let the statistic Sl^\p,u,e) be as in (l2TI) . cor- 

i'l') (2) 

responding to f^ J. Then, in the limit t the process S, T/i, n, £), n G (0,1), converges in law 
to the centered gaussian field S^^\p,u,e) having the asymptotic covariance 


( 2 )/ 


Cov 




_ j-/i(31og£-K: + log|M-v|), if\u-v\ >£, 
|-/i(41og£-2fc), ifu = v, 

+ 0{e). (148) 


Proof This follows from Theorem [24] and the calculation of the scalar product in Lemma [431 | 


On the other hand, it is elementary to see from (iTTll -dT^ that the gaussian free field at scale £ 
satisfies 

—plog\u — v\, if |m —v| > £, 

/i(l — log£), if M = V, 


Co\(^(0^^e{u),(0^,e{v)j = I 

so that the centered gaussian free field at scale £, 

= tU^_e(M) — ft)^^g(O), 

satisfies 

Cov(^d)^,e(M),d)^,e(v)^ = 


(149) 


—/i(log£ — 1 +log|M — v| — loglnl — log |v|), if |m — v| > £, 


— 2/1 (log £ — 1 — loglnl), if M = 


(150) 


(151) 
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As we remarked in our discussion of (l95l) . the choice of truncation, and, in particular, the choice of 
the constant, e—independent term in the covariance of the free field, has no effect on the law of the 
total mass so long as the key normalization condition in (IMl) holds. Hence, for our purposes, we can 
re-define g(M) to be 


Cov [( 0 ^^e{u), 0)^,6 (v)] 


{ —/I login — v|, e < |m — v| < 1, 
At(fc-log£), 


(152) 

(153) 


so that ( 0 ^^g{u) has asymptotically the same covariance as that of in (I145I ). Hence, 




inlaw 

~ 0)^,6 ( n ) 


(154) 


as stochastic processes in the limit e —)• 0, up to zero-mean corrections having covariance of the order 
0{e). Similarly, by comparing covariances of (ft,£) in (11481) and (0^^e{u) in (11531) and recalling 
that 5(2)(jU,M,£) is centered, while the mean of ( 0 ^^e{u) is given in (11521) . we obtain, also up to 0{e), 

5(2)(^,m,£) ‘~"'cu^ g(M) + ^(^-^(log£-K-),-;U(31og£-fC)), (155) 

where c/K(—(/i/2)(log£ — fc), —/i(31og£ — fc)) is an independent gaussian random variable having 
the mean — (/r/2)(log£ — k) and variance —/i(31og£ — k). 

We can now explain the origin of Conjecture 12.31 The basic idea is to form the exponential 
functional of the statistic s\'\lJL,u,e) and compute its Mellin transform by analogy with the gaussian 
free field so as fo obfain fhe fofal mass of fhe limif lognormal measure in fhe limif. The principal 
obsfacle in fhe firsl case is fhaf behaves like d)^ g(M) as opposed fo (0^^e{u) so fhaf ifs 

exponential functional does nol exisf in fhe limif £ —)• 0. This obsfacle is overcome by appropriately 
rescaling fhe Mellin Iransform as shown in fhe following fheorem. To Ibis end, we will firsl formulaic 
a general proposition and Ihen specialize if fo Coniecfure 12.31 In whal follows (p{u) can be a general 
lesl funclion, however for clarify, we will reslricf ourselves fo 

(p{u) = u^'{l-uf^,XuX2> (156) 


Theorem 4.5 (Rescaled Mellin transforms) Let 0 < < 2 and 1 be a subinterval of the unit inter¬ 

val. Then, for Re{q) <2jn we have 


For —(1 + Ai)//i — 1/2 < Re(^) Kljllwe have 


= E 


lim e 

£->0 


,M(logg-K-)2(^ 


E 




= E 


u)M^{du 


(p{u)Mn{du 


Moreover, the positive integral moments of order n <2/p satisfy 

\ n(n+l) 


e->0 


lin^gM(ioge->c)—E / = m[(p(5;)]]^|5,--5^-|-^'r/5i...r/5„, 

I '=1 '■<; 


e->0 


i=l 


KJ 


(157) 


(158) 


(159) 

(160) 
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The proof requires two auxiliary calculations. The first involves a version of Girsanov-type theorem 
for gaussian fields, which exfends whaf we used in our original derivafion of infermiffency differenli- 
afion, cf. Secfion 8 of (431. For concreteness, we sfafe if here for fhe G)^,e(f) process. 

Lemma 4.6 (Girsanov) Let 1 be a subinterval of the unit interval. Then, we have for Re(p) 




e^O 


= E 


u)M^(du 


(161) 


andforRe{q) > —(1 + Ai)//X — 1/2, Re(p) < min{2//i, 2//r(l + Ai) + 1 +2Re(^)}, 


lim e 

e->0 


ii{\oge-K)Slspl 


—E 




= E 


u^‘^ q){u)M^{du 


(162) 


Proof The integral JiU^‘^ (p{u)M^{du) is defined for jj. Re(^) + X\+ pt /2> —\ and ifs Mellin frans- 
form for Re(p) < min{2//i, 2//r(l + Ai) + 1 +2Re(^)}, cf. (11371 ). Nexf, we will esfablish fhe 
idenfifies 




~E 


-qWfifO) 


-M 


(p{u)e 




du 


g-qoigA^)l I 


= E 
= E 






(163) 

(164) 


This resulf is equivalenf fo a change of measure. Infroduce an equivalenf probability measure 

d.S = e-M{'f-loge)^g-9%,£(0) 


(165) 


where is fhe original probabilify measure corresponding fo E. In facl, if follows from (11521 ) and 
(11531) fhaf 


E 


^-AK-\oge)Sl^^-qmnAO) 


= 1 . 


(166) 


Then, (he law of (he process u —)■ (o^^e{u) — qCoy[(0^^e{u), ft)^,e(0)) wifh respecf fo tlP equals (he 
law of (he original process u —)• ( 0 ^^e{u) wifh respecf fo .S. Indeed, if is easy fo show fhaf fhe fwo 
processes have fhe same finife-dimensional disfribufions by computing (heir characferisfic funcfions. 
The confinuify of sample pafhs can (hen be used fo conclude fhaf fhe equalify of all finife-dimensional 
disfribufions implies (he equalify in law. Once (his equalify is esfablished, (hen 




= E 
= E 








(167) 


which completes (he proof of (11631) and of (1161b by faking fhe limif e —)• 0. The argumenf for (1162b 
is essenfially (he same excepf fhaf one needs fo resfricf p so fhaf fhe righf-hand side is well-defined. 
Given (1117b . (he conditions are Re(p) <2l\l and 1 — Re(p) -|-2Re(^) + { \ + Xf)!/pt > 0. | 


The second step in fhe proof of Theorem 14.51 enfails a key momenf calculation, in which we 
compute (he asymptotic of fhe rescaled positive infegral momenfs in terms of generalized Selberg 
integrals. 
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Lemma 4.7 (Single moments) Let n <1l\i and 1 be a subinterval of the unit interval. Then, as 

£ —)• 0 , 


E 




|5(—5y| ^ds\---dsn, (168) 


(=1 i<i 


E 


(p(M)e“f'’'=(“) du 


e 


/i(-loge+K') 


n(n+l) 


ri['p(^'>r]n |5, —iyl ^dsi-'-dsn- (169) 


/« 


!=1 


i<] 


Proof We will give the proof of (I169I ). The idea is to discretize as in (|95] ). Let N = |/|/£, COj = 
d)^,e{^j), ‘5'; = £7, j = 1 • • W. As d)i, • • • , d)A? are jointly gaussian with zero mean, we have 


N 


N 


i=l jl-jn = t 


(170) 


Now, using (11521) and (11531) . 

n n n 

Var(d)y, +... + d)yJ = £ Vard)y,+2 £ Cov(d)y,,d)y,) =2/i£(-log£ + fC + log5yJ + 


k=\ 


k<l=\ 


k=\ 


+ £ (-loge + K:-log|5y,-5y,|+log|5yJ+log|5y,l), 


k<I=\ 


= 7t(-log£ + fC)n(n+l)+2/In£log5y^-2/r £ log |5y^ - i'y, |. (171) 

k=\ k<l=l 


It follows from (11701) that 

E[£f £ n[,>(»,-j^"] n (172) 

i=\ y,...y„=U=l k<l=\ 

In the limit we have —>• fidsk, hence the result. The proof of (11681) follows by re-labeling 

q){u) ^ u^^"q){u). I 

We note that the same type of argument gives a formula for the joint moments. 

Lemma 4.8 (Joint moments) Let N = n + m, N <2/jx, and f f 2 be subintervals of [0,1]. Then, as 

£ —)• 0 , 


E 


(^ j (^ j (p 2 {u) duj 

h h 




M(-loge+)c)M 


/"x/^ 


fi[<Pi(^c)^n 


X n ^dsi---dsN. 

i=n+l i< j 

(173) 


We can now give a proof of Theorem 14. 5 1 and explain the origin of Conjecture 12.31 
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Proof of Theorem I4.5[ Let Reip) < 2/p. By Lemma lL6l 


£->0 


-gcon,e(0)f 


= E 


(p{u)M^{du 


(174) 


In particular, (11741) holds for ^ = p, if Re(^) <2jp. On the other hand, the left-hand side of (11741) 
can be reduced when q = p as follows. Substituting the definition of a)^^e{u), cf. (11501) . we obtain 


lim^p(logr-c)^E 

£->0 





g-\ 


e^O 



(p{u) 


(175) 


The argument for (11581) is essentially the same except that we use (11621) instead of (11741) . Given 
Re(^) > —(1 + Ai)//i — 1/2, Re(p) < min{2//i, 2jp{\ + Ai) + 1 +2Re(g)}, if <7 = p the condition 
Re(p) < min{2/jU, 2//i(l + Ai) + 1 +2Re(^)} is equivalent to Re(g) > —(1 + Ai)2/jU — 1 and is so 
satisfied. | 


Proof of Conjecture 12.31 Lef (p{u) = 1 and I = [0, 1] in Theorem 14.51 Recalling fhe asympfofic 
equably of laws of d)^^e{u) and Sf^\p,u,e) al finite e > 0 in fhe limil f —)■ 00 thaf we esfablished 
in (11541) above, we formally inlerchange (he order of t limil and u infegrafion, lake (he t limit by 
Corollary 14.21 then replace S^^\p,u,e) with d)^^e{u) by (11541) . and, finally take the limit e —)• 0 by 
Theorem 14.51 The right-hand sides of these equations are known by Conjecture I3.6l and described in 
(11171) . with Al = A2 = 0, and Ai = pq, A2 = 0, respectively, cf. Theorem 13.21 above, hence (l49l ) and 
(l52l) . For example, the argument for (|4^ is as follows. 



limE 


£->0 

yoo 

L \Jo J J J 



E 


£->0 


LVjo '' J 


= lime(^(i°g"-'^) 
£->0 


g(g+l) 

2 


E 





= E 



(176) 


( 2 ) 

The derivation of the weak conjecture for S) (p,u,e) follows the same steps but is simpler as 

f2') 

it is based directly on the asymptotic equality of laws of co^^e{u) and S) {p,u,e) at finite £ > 0 in 
the limit t ^ oo^ cf. (11551) above. By taking the exponential functional of 5) {p,u,e) as in (|4^ and 
using (11551) . we obtain the Mellin transform of the total mass of the limit lognormal measure times 
the Mellin transform of exp(M^(—(/i/2)(log£ — fc), —/r(31og£ — k))). The latter accounts for the 
scaling factor in (|4^ . and the result then follows from Conjecture 13.61 with Ai = A2 = 0. Similarly, 
(|5TI ) follows in the same way from Conjecture 13.61 with Ai = pq, A2 = 0. | 


We next proceed to the derivation of Conjecture 12.41 The basic idea is to let e{t) approach zero 
“slowly” compared to the growth of t to infinity so that the statistic {p, u) in (l3^ behaves as the 
centered gaussian free field for / = 1 or the gaussian free field plus an independent gaussian for i = 2 
at the scale e{t). The main technical challenge of quantifying the required rate of decay of e{t) is 
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that neither Theorem 12. 1 1 nor IT2] applies to the statistic Theorem 12. II does not apply as the 

variance of {l^,u) is divergent in the limit f —)• oo. Theorem l2.2l does not apply as the asymptotic of 
the variance \w\\Xu{w)\^ dw oc log (log t/A (t)) in (jT]) is different from the asymptotic of our 

oo __ 

variance / \w\\{Xu* <j>e){w)\'^ dw o<^ — loge, cf. Lemmas |4. 1 1 and 1431 above. Instead, we need a slight 

— oo 

modification of Theorem 12.21 Lemma WM shows that for a sufficiently slowly decaying £{t) one can 
obtain a limiting gaussian field having e(f)-dependent asymptotic covariance. Lemma |4. 1 01 shows 
further that the limiting covariance can be approximated by the scalar product in Theorem 12.11 Let 
St{f) be as in (O. For simplicity, we restrict ourselves here to finite linear combinations of indicator 
functions Y.CkXut, 0<Uk<l. 

Lemma 4.9 (Modified convergence) Let X{t) satisfy (EJl and ef) satisfy (1^ for i = I or (1^ for 
i = 2, respectively. Let Xu\x) = X[o,u\{x) and x^\x) = A[-i/e(0,«] W- Define 


\ogt/X(t) 




w 




dw, 


-log?/A (0 


(177) 


then, as t ^ 

(178) 

The proof is sketched in the appendix. Let fel{x) = {Xu'^ * ^e){x) be as in Lemmas 1411 and 1431 It 
follows by linearity that the covariance has the asymptotic form as t —)• oo 

logf/A(7) _ _ 

c<-v( 4 (/i;!,,.). 5 ,(/i;!,,))=Re / (179) 

-logr/A(() 

It remains to show that this asymptotic is the same as that of the covariance in (11391) and (11401) for 
i = 1 and (11461) and (11471) for i = 2. This is shown in the next lemma. 

Lemma 4.10 (Limit covariance) Let X{t) satisfy (Ill) and £{t) satisfy (l3^ for i= \ or (1351) for i = 2, 
respectively, and the scalar product be as in (ID). In the limit t —)• 00 ^ 


log7/A(() ^^ -- 

-log7/A(7) 

Proof of Lemma 14.101 Given the equality between the representations of the scalar product in 
dH) and dSl), we need to estimate Jw)/^||j^^(>v)(iw. Denote a = \ogt/X{t), 

fe{x) = fei{x), and ge(x) = felfx). Then, by substituting the definition of the Fourier transform and 
integrating by parts. 


Re J \w\fe(w)ge(w)dw = 


f'e {^) s'e (y) (cos (rvy) - cos (w (y - x))) dxdy 


(181) 
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The cos(>vj) term can be replaced with cos(>v) when integrating from 0 to a and dropped when 
integrating from a to infinity because f f'(x)dx = 0. Using the Frullani integral in the form 


oo 



0 


-cos(w(y-x))) =log|y-;c|, 


(182) 


we wish to show that the remainder term satisfies the estimate in the limit a —)■ oo and £ —)■ 0, 


oo 




a 


/e W §^ 3 ") cos (w(y - x) )dxdy 


0 (l/a£). 


(183) 


Let i = 1 for concreteness. Recalling /e(x) = <pe{x) — <pe{x — u), gg(y) = 0e(j) — <Pe{y — v), and the 
assumption that (p (x) is compactly supported so that the only values of x and y that contribute to these 
integrals are concentrated in intervals of size £, it is sufficient to show 


J ~ J peix) cos (wx)dx = J— J p{x) cos (wx)d. 


ae 

=0{\/ae) 


(184) 


which follows by integrating the inner integral by parts. | 

Remark One concludes that the statistic s|^^(/x,m) and the centered gaussian free field d)^g(f)(M) 
have asymptotically the same law, as do the statistic S] {li,u) and the gaussian free field plus an 
independent gaussian, as in (11541) and (11551) with the scale e{t). The difference between 
and Sf {ll^u) on the one hand and the corresponding free fields on the other is that these statistics 
are only gaussian in the limit f —>• oo^ whereas the free fields are gaussian for £ > 0. We will assume 
that deviations of these statistics from their gaussian limits can be ignored for sufficiently large t for 
the purpose of carrying out calculations. We can quantify this assumption in terms of the behavior of 
near its gaussian limit in (11781 ). by conjecturing the asymptotic of St‘\iJ.,u) to be 

Sl'\ix,u) ' ='" a/'V(0,1) +o(l), (185) 

which means that converges to ^{0, 1) faster (in the sense of variance) than oc — log£(t) 
goes to infinity. We can heuristically estimate the rate of convergence of f/'^ to ^(0, 1) by the 
relative error in the formula for the variance of our statistics. It follows from (11831) in the proof of 
Lemma |4.101 that the relative error is 0(A(t)/logt£(t))/|log£(t)| so that the rate of convergence 
^1/1 log £ (f) I by (l34l) as required for (11851) to be true. Obviously, the assumption in (11851) is a major 
gap between the weak and strong conjectures. 


Proof of Conjecture 12.41 The argument is the same as that for Conjecture 12.31 above except that 
we use Lemma 1431 instead of Theorem 12. II and Lemma |4. 101 instead of Lemmas |4. 1 1 and 1431 and, in 
addition, assume that (ft,«) is near-gaussian for large but finite t as remarked above in (11851) . | 


The proofs of the corollaries are given below. They are the same for either the single or double 
limit on the left-hand side. For concreteness, they are stated for the single limit as in Section 2. 
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Proof of Corollary 12.51 We computed the positive moments and verified (l57l) in Lemma 1471 On the 
other hand, using (l49l) . (l5^ follows from the known formula for the positive integral moments of the 
Selberg integral distribution, cf. (I118I) . The argument for (|5^ and (1^ is the same. | 

Proof of Corollary 12.61 This again follows from Conjecture 12.41 and the known formula for the 
negative integral moments of the Selberg integral distribution, cf. (II 191 ). | 


Proof of Corollary 12.71 This is a special case of Lemma 14. 8 1 | 


Proof of Corollary 12.81 (1^ is a special case of (11221) corresponding to Ai = A 2 = 0 . 
case of the following asymptotic expansion in the limit /r —)■ 0 , 


is a special 


931 (^ 1 /i, 7 tz, 0 ) ~exp(£(^^^ 


r=0 


|- fir+2(^+1)+^r+2(^) ^ 


+ 


Br+liq - 2 z) - Br+2i-2z) - 26^+2 
r + 2 


r+ 1 r + 2 

Br+2 {q - I- 2 z)- Br +2 { 2 q-\- 2 z)- 


-^] + (C(r+l)-l)[ 


r + 2 


(186) 


This expansion coincides with the intermittency expansion of E 




as one can verify 


by the methods of the previous section. Another way of deriving (11861) is to follow the argument that 
we used in Il48l to derive (11221) from (11171) . Then, (l 6 ^ corresponds to z = <7 in (11861) . | 


Proof of Corollary 12.91 We will give the argument for the second statistic. We have by (11551) . 

ii+A iS2+A ii+A i2+A 

Cov^log J log J «Cov(log J log J 

.?! ^2 ^1 ^2 

-At(31og£(0-fc), (187) 

and the result follows by (IMt and the stationarity property of the limit lognormal measure. The 
argument for the first statistic is similar but is more involved as it requires a generalization of (IMl) for 
the centered gaussian free field, which follows from (11731) in the same way as ([98]) follows from (1971) . 
The details are straightforward and will be omitted. | 


Proof of Corollary 12.101 Given Conjecture 12.41 these results follow from Theorems 13.21 and 13.41 

I 


Proof of Corollary 12.111 Given the multifractal law of the limit lognormal measure, cf. (|9T]) . it is 
sufficient to show the identity for any 0 < a < 1 


limeM(logr(0-r)2(^E 

t-+oo 


-M, 




= E 



(188) 


which is a special case of (11571) corresponding to (p{u) = 1 and 1 = [0,a] (by formally replacing 
d)^,e(w) with (ft,«) as in the derivation of Conjecture 12.41 abovel. | 
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Proof of Corollary 12.121 Given the multiscaling law of the limit lognormal measure, ^ follows 
from (l92l) . To prove (1751) . we need to generalize (l9T]) to the following identity in law, 


Jo Jo 


(189) 


where Q.s is as in (I86b and (1871) . which implies the multiscaling law 


E 


u^‘^M^{du 




E 


oc 


(190) 


as a function of 5 < 1. Finally, (11891) is a simple corollary of 


Remark It should be clear from Theorem 13.21 and Theorem 14.51 that one can obtain general explicit 
formulas for q){u) = u^'{ \ — u)^-. The reason that we restricted ourselves to Ai = A 2 = 0 in Conjec¬ 
tures 12.31 and 12.41 is simplicity. The reader who is interested in the general case can easily find the 
desired formulas in Section 3. 


5 Conclusions 

We have formulated two versions of a precise conjecture on limits of rescaled Mellin-type transforms 
of the exponential functional of the Bourgade-Kuan-Rodgers statistic in the mesoscopic regime. The 
conjecture is based on our construction of particular Bourgade-Kuan-Rodgers statistics of Riemann 
zeroes that converge to modifications of the centered gaussian or gaussian free fields. The sfafisfics 
are defined by smoofhing fhe indicator funcfion of cerfain bounded or unbounded subinfervals of fhe 
real line. The smoofhing is effecfed by a rescaled bump function. In fhe weak version of fhe conjee- 
lure fhe asymptotic scale of fhe bump funcfion e is fixed so lhaf fhe resulting slalisfics 
i = 1,2, satisfy fhe Bourgade-Kuan-Rodgers Iheorem. In fhe slrong version, fhis scale e{l) is chosen 
to be mesoscopic, X{t)/\ogt <C £(t) 1 for i = 1 and bolh A(t)/logt <C e(t) <C 1 and E{t) S> 1/A(t) 

for i = 2, so lhaf fhe slalisfics sl’\lJ,,u) satisfy an extension of fhe Bourgade-Kuan-Rodgers Iheorem 
thal we formulafed in fhe paper. We have compufed fhe double limil over fhe scale e and t in fhe weak 
case and fhe single limil over t in fhe slrong case of rescaled Mellin-lype Iransforms of fhe exponen- 
fial functional of bolh slalisfics as if fhe slalisfics were fhe cenlered gaussian free field or fhe gaussian 
free field plus an independenf gaussian random variable, respectively. The exponential functional of 
fhe gaussian free field is an imporlanf objecl in mafhemafical physics known as fhe limil lognormal 
stochastic measure (or lognormal mulliplicafive chaos). By using a Girsanov-fype resulf, we have 
found an appropriafe rescaling faclor to compute Iwo Mellin-lype Iransforms of fhe exponential func- 
lional of fhe cenlered field in lerms of fhe Mellin Iransform of fhe exponenlial functional of fhe free 
field ilself, i.e. fhe Mellin Iransform of fhe lolal mass of fhe limil lognormal measure, resulling in 
fhe conjeclure for fhe firsl slalisfic. The conjeclure for fhe second sfafislic follows direcfly from ifs 
convergence to fhe gaussian free field plus an independenf gaussian, fhe laller being responsible for 
rescaling. In bolh cases, fhe rescaling factors are determined by fhe asymptotic scale and fhe choice 
of fhe bump funcfion lhaf effecl fhe smoofhing. Finally, our conjeclural knowledge of fhe disfribufion 
of fhe lolal mass of fhe measure has allowed us lo calculale a number of quantifies lhaf are associaled 
wifh fhe slalisfics exaclly. 


31 










The principal difference between the weak and strong conjectures is their informational content. 
Conjecture 12.31 alone can be thought of as a number theoretic re-formulation of Conjecture 13.61 on 
the equality of the Selberg integral distribution and the law of the total mass of the limit lognormal 
measure. It associates the Selberg integral distribution with the zeroes but does not contain any infor¬ 
mation about their distribution that is not already contained in the Bourgade-Kuan-Rodgers theorem. 
Indeed, the order of the u integral and the t limit can be interchanged at a finite £ > 0, resulting in the 
Mellin transform of the exponential functional of a gaussian field, which converges fo fhe centered 
gaussian free field when / = 1 or fhe sum of fhe gaussian free field plus an independenf gaussian 
random variable when / = 2 as £ —)• 0, so fhaf fhe weak conjecfure only requires fhe t —)• oo limif of 
fhe S^’'^ (/I, M, £) sfafisfic. On fhe ofher hand, Conjecfure [2]4] combines fhe £ —)• 0 and t —)• oo limifs info 
a single limif so fhaf fhe order of fhe u infegral and fhe resulting limif can no longer be inferchanged 
because fhe slafislic unlike s\'\lJL,u,e), converges fo fhe centered gaussian free field when 

/ = 1 or fhe gaussian free field plus an independenf gaussian when / = 2 af our mesoscopic scale 
E{t) and so becomes singular as f —oo. The sfrong conjecfure assumes fhaf deviations of s\^\}X,u) 
from ifs gaussian limif are negligible af large buf finile t, fhereby providing some new informalion 
abouf fhe sfafisfical disfribufion of fhe zeroes af finite t. In particular, as fhe sfrong conjecfure fifs info 
fhe framework of mod-gaussian convergence, fhe normalify zone and precise deviations of fails of 
our exponential functionals can be compufed from fhe general fheory and explicif knowledge of our 
limifing functions. 

We have provided a self-conlained review of some of fhe key properties of fhe limif lognormal 
measure and fhe disfribufion of ifs fofal mass fo make our work accessible fo a wider audience. In 
particular, we have covered fhe invariances of fhe gaussian free field, fhe mulfifracfal law of fhe limit 
measure, the derivation of the law of its total mass by exact renormalization, and several character¬ 
izations of the Selberg integral distribution, which is believed to describe the law of the total mass. 
The Selberg integral distribution is a highly non-trivial, log-infinitely divisible probability distribution 
having the property that its positive integral moments are given by the Selberg integral of the same 
dimension as the order of the moment. We have reviewed both its analytic and probabilistic structures 
that are relevant to out calculations. 

We have provided a number of calculations that support our conjecture. Our calculations are 
universal as they apply to any asymptotically gaussian linear statistic having the covariance structure 
that is given by the Bourgade-Kuan-Rodgers formula. In particular, they apply to the GUE statistics of 
Fyodorov et. al. Il27]l and provide a theoretical explanation, via convergence to the gaussian free field, 
for why our resulfs for fhe Bourgade-Kuan-Rodgers and GUE sfafisfics are fhe same. Our argumenfs 
are however nol mafhemafically rigorous as we do nof know how our sfafisfics behave near fheir 
gaussian limifs. Our assumpfion fhaf fheir behavior in fhe limif can be used fo do calculafions near 
the limit in the strong case is the principal mathematical gap between the weak and strong conjectures 
that renders our use of the free fields in place of fhe sfalislics heuristic. We have quanlified fhaf fhe 
variances of our sfafislics should converge fo fheir asymptotic limits faster than l/|log£(t)| for this 
assumption to be valid and explained heuristically why we expect this to be true. 

In broad terms, on the one hand, our conjecture relates a limit of a statistic of Riemann zeroes 
with the Selberg integral and, more generally, the Selberg integral probability distribution and so 
associates a non-trivial, log-infinitely divisible distribution with the zeroes. On the other hand, our 
conjecture implies that the limit lognormal measure can be modeled in terms of the zeroes. As this 
measure appears naturally in various contexts that involve multifractality, we can speculate that there 
is a number theoretic interpretation of multifractal phenomena. In particular, a proof of (even the 


32 


weak) the conjecture might lead to a number theoretic proof of the conjectured equality of the law of 
the total mass of the limit lognormal measure and the Selberg integral distribution. 

We have interpreted our statistics as fluctuations of the smoothed error term in the zero counting 
function so that our conjecture gives rescaled Mellin transforms of the exponential functional of these 
fluctuations. Aside from verification or disproof of our conjecture, it would be quite interesting to see 
what properties of the error term follow from the conjecture and to compute the rate of convergence 
of our statistics to their gaussian limits. We believe that our conjecture is only valid for 1 <C A(f) <C 
logf and breaks down for A(t) ~ 1 due to the expected presence of yet-to-be-determined arithmetic 
corrections. Finally, the appearance of the Selberg integral distribution in the work of Fyodorov and 
Keating OOll and in our paper begs the question of formulating a general statement about logarithmic 
correlations in the statistical value distribution of log which we believe will clarify the relationship 
between their conjecture and ours. 
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A Appendix 

In this section we will give a proo^ of (IThl) and a sketch of the proof of Eemma 14.91 The starting 
point is the second equation following Eq. (16) in ifTSl . 




1 


\Y{t)-{x + iy) 

Recalling the identity 


Re 




1 log 1 1 

2 X^^X^)"“ C ^2 ' X{t) 




a{ty 

(A.l) 


lim 

>>->0 


x)Im 


dx = 7if{Y), 


^Y-{x + iy) 

we multiply (lA.ll) by f{x) jn, integrate with respect to x, and let y —)• 0. We then obtain by Q 


(A.2) 


St{f) 



dx+ 0(l/A(i)), 
dx + 0{\/X{t)y 


= - J fix) s(^(ot + dx + Oil/Xit)). 

R ' 


(A.3) 

(A.4) 

(A.5) 


It should be noted that Re(^'/^)(5) is integrable along the critical line due to Hadamard’s factoriza¬ 
tion, which shows that the singular part of /Qf) near the zeroes is of the form \/{s — y) and 

^*This proof is due to Nickolas Simm. We originally established this result as a corollary of Eq. (17) in ca Upon 
seeing our derivation, Nickolas Simm found a much simpler proof, which is given here. 
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so is purely imaginary along the critical line, hence vanishing upon taking the real part. It is clear 
from (lA.ll) that the error term is indeed C?(l/A(t)) provided f{x) is compactly supported and the 
support is independent of t. If the support depends on t as in (1^ . provided its size<C tX{t), the error 
is 0[\\f\\i/X{t)) as indicated in footnote 7. 

The proof of Lemma 149] follows verbatim the proof of Theorem I2.2I given in lITSl . We will only 
indicate the main steps here and refer the reader to ifTSl for all details, including definitions of special 
number theoretic functions that are needed in the proof. 

Proof of Lemma 14.91 Let Xu{x) = X[o.u]{^)j * = 1- Let {cu} and {u^} be fixed and denofe 

/fW - (A.6) 

Then, ft and // are bounded in uniformly in t, and //' satisfies 

||/;'||i=O(l/£(0). (A.7) 

The Fourier fransform of f satisfies uniformly in f as w —)• oo 

w|7f(w)p, (w|7f(w)p)' = 0(l/w). (A.8) 

Then, by fhe approximafe explicif formula of Bourgade and Kuan, cf. Proposition 3 in lITSl . 


Stift) = 



(A.9) 


where A„(n) denofes Selberg’s smoofhed von Mangold! function. In our case, as in fhe case of 
Theorem l2.21 fhe error ferm is of fhe order 0{X{t)/\ogt e{t)). Lef Ot be as in (11771) and denote 


A V2 y A^(p)^/logp 



Then, by (IA.9! ) and Lemma 4 of ifTSl . 


(A. 10) 


Stif) 



+ 0 


m \ 

logte(t)/ 


+ «( 1 ) 


(A. 11) 


in fhe limif t —)• oo. Finally, by Lemma 5 and Proposifion 6 of ifTSl . Yf 
complex normal variable 


Y,^ 


M^(0,1) + /aL(0,1) 

7! 


converges fo fhe sfandard 
(A. 12) 


The resulf for i = 1 follows. 

Now, lef Xu{x) = X[-i/e{t),u] /, be as in (IA.6I) . The argumenf given above still goes fhrough 

buf requires somewhaf more delicate esfimafes. We have for n = 0,1,2, 


||/W||i = o(i/^^(0), 

|xlogx//"^||i =0{log{\/e{t))/e^), 


(A. 13) 
(A. 14) 
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so that (IA.9I ) still holds with the error term of the order 0{X{t)/logt e{t)). To verify (lA.llI) we need 
Lemma4 of ifTSl . which requires that ||/j||i be uniformly bounded, whereas we have (IA.13I) instead. 
However, a careful reading of the proof of Lemma 4 indicates that it still holds provided 


which in our case translates into 

1 

, , . . , = o(l). 
e{t)X{t) 

Thus, (lA.l 111 holds given the condition in (l35l) . Similarly, 

(IA.8I) . whereas we have instead 

w\l{w)\^ = 0{l/w), (A.17) 

{w\l{w)\^)' = 0{l/we{t)). (A.18) 

Once again, a careful reading of the proof of Lemma 5 indicates that the bound in (IA.18I) is sufficient 
provided \/e{t)X{t) = o( 1) as in (l35l) . The rest of the argument goes though verbatim. | 
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